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ABSTRACT 


Pressure fluctuations due to the interaction of a shock wave with a turbulent boundary layer 
are investigated. A simple model is proposed in which the shock wave is convected from its 
mean position by velocity fluctuations in the turbulent boundary layer. Displacement of the 
shock is assumed limited by a linear restoring mechanism. Predictions of peak root mean 
square pressure fluctuation and spectral density are in excellent agreement with available 
experimental data. 
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1.0 INTRODUCTION 


An important- source of surface pressure fluctuations on high speed aerodynamic vehicles Is the 
oscillation of shock waves. Figure 1 shows the static and fluctuating pressure levels associated 
with supersonic flow ahead of a flare or a three-dimensional protuberance. Other basic flow 
geometries containing oscillating shock waves exist, such as oscillation of a near normal 
terminal shock in transonic flow. Reference 1 provides detailed descriptions of fluctuating 
flow fields, and contains a comprehensive review of available experimental data. The flow 
field of Figure 1 will be briefly described here to provide a framework for the shock oscillation 
model presented in the next section. 

There are three basic sources of fluctuating pressure in the flow field shown in Figure 1: the 
attached turbulent boundary layer ahead of the shock; the separated region behind the shock; 
and the oscillation of the shock itself. Figure 2 shows spectra of pressure fluctuations in these 
three parts of the flow field ahead of a 45° wedge, from Reference 2. Robertson (Reference 1 ) 
has pointed out that the spectrum at the shock location represents the spectrum associated with 
an oscillating shock wave alone, plus some fraction of the attached and detached spectra, as 
these two environments alternately appear at the mean shock location as the shock oscillates. 
Fluctuations from the attached boundary layer are much smaller than from separated flow and 
may be neglected, so the spectrum of fluctuations at the mean shock location is given by shock 
oscillation and separated flow environments. In Figure 2, the spectral shape follows the separated 
flow spectrum for f 6 / u^ > 0.2; the spectrum for shock oscillation alone would follow the 
dashed extrapolation. 

Much of the analytical work on oscillating shock waves has followed the approach used by 
Trilling (Reference 3) for interaction of a shock wave with a laminar boundary layer. In that 
analysis, a harmonic perturbation to the mean flow was assumed, and it was shown that oscilla- 
tion could be self-sustaining for certain frequencies. This approach is felt to be unsatisfactory 
to the present problem for two reasons. First, the source of the initial disturbance is not 
identified. In the experimental investigation of Reference 4, a spark discharge was used to 
artificially stimulate oscillation. Second, Trilling's analysis would lead us to expect dominant 
frequencies in the spectrum. The spectrum of Figure 2 does not contain any such peaks. 
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Figure 1 . Composite Schematic of Protuberance Flow Field Characteristics, 

M = 1.60, h/D = 2.0, from Reference 1 
oo 
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O Attached turbulent boundary layer 

^7 A Non-homogeneous separated flow 

□ Shock wave oscillation 

Extrapolated shock oscillation alone 

Figure 2. Longitudinal Distribution of Pressure Fluctuations and Typical 
Power Spectra in Vicinity of Supersonic Flow Separation 
Ahead of a 45° Wedge, from Reference 2 
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The broad band nature of observed shock oscillation spectra in the presence of a turbulent 
boundary layer suggests the following model. The mean location of the shock is a stable 
position, governed by mean flow conditions. (If this were not so, the shock would not be 
there.) As each turbulent eddy passes, the shock is convected upstream or downstream. 

The stability of the mean shock location causes a limit to the excursion distance. In the next 
section, an equation of motion for shock location is postulated, based on convective displace- 
ment and a linear restoring mechanism. Statistical properties of the shock motion are found 
from which mean square pressure fluctuations and spectra are calculated. 
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2.0 ANALYSIS 


2.1 Equation of Motion of the Shock Wave 

Velocity fluctuations in the x-direction (parallel to the wall) at a fixed point within the 
boundary layer may be represented as 

u T = a € p (t) (1 ) 

I 00 

X 

where p( t) is a random funcHon of time with < p 2 y = 1 , and € is the turbulent Mach 

l 

number < u 2 } 2 /a^ # ^ one-dimensional model is adopted, neglecting velocity fluctua- 
x 

tions in other directions. The shock wave within the boundary layer is then convected in the 

x-direction with speed u^ . 

x 

After the shock wave is displaced a distance x (x small compared to scale length of the flow 
geometry), the flow geometry will be disturbed by an amount proportional to x . This in turn 
will disturb the pressure field by an amount proportional to x . The shock wave will then move 
with velocity proportional to x . Since the shock will move back toward its original position, 
this restoring velocity is given by 


u . = - j3 x (2) 

restoring 

where (3 is a constant depending on the flow geometry. This constant will be deduced as the 
analysis proceeds. 

The net speed of the shock wave is the sum of Equations (1 ) and (2): 

u = a £ p (t) - !3 x (3) 

CO 

The equation for the shock location is 
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(4) 


a? + ' Jx = e ° ro tJ <f) 

Equation (4) may be integrated to give 



(5) 


In the following subsections, various statistical properties of x are calculated, leading to 
predictions of the intensity and spectral distribution of pressure fluctuations. 


2J .1 Mean Square Shock Displacement — Squaring Equation (5) and taking the ensemble 
average. 


<x 2 > = € 2 a 2 e- 20 
00 


’[[ <H(£ )h(£ )>e 3(5 l +5 2> d5 t (6) 

/o ^0 1 


Changing integration variables from | | to center of mass and relative coordinates 

+ £ )/2 , £ “ £ ( / it is straightforward to obtain the following equation for t 55> 

v I 2 2 \ ' 


T : 


<X 2 > = € 2 a 2 (1 - e -' 3t ) 

00 p 


(7) 


where: 


T = 




R 

H 


M 2 > 

< M 2 > 


If is assumed that x = 0 af f = 0. 

The calculation leading to Equation (7) is identical to the one given in Reference 5 for the 
Brownian motion of a particle with linear damping. Details of the calculation may be found 
in the next section, where the correlation function is obtained. 
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( 8 ) 


For large time, t » -j- , <x 2 > reaches a finite limit: 


-<^> . «... i 


Because of this constant limit, x for large time may be treated as a stationary random function. 


It is interesting to compare Equation (8) with the well known result for 0 = 0 , corresponding to 

an unrestricted random walk: 

^ 



2 € 2 a 2 T t 

00 p 


(9) 


Without the restoring mechanism, displacement continues to grow with time. The asymptotic 
value given by Equation (8) is also given by Equation (9) when t = 1/2/3 . If the restoring 
mechanism is thought of as limiting the time during which an unrestricted random walk takes 
place, then it is clear that 0 must be related to the integral time scale of the x motion. The 
integral time scale is found from the correlation function, which will now be calculated. 


2.1 .2 Correlation Function and Integral Scale — The correlation function (t) is defined 
as: 


R 

x 


(t , t ) 


2 


<x(t ) X (t ) > 
1 2 

<x 2 > 


00 ) 


x is assumed to be homogeneous and isotropic, so that 


R (t , t ) = R (t -t ) = R (t -t 
x l 2 x i 2 x 2 1 


Mr) 


OD 
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Working from Equation (5) and setting < x(t ) x(t ) > = <xx> 

1 2 


y x -fit -fit 2 
<xx> = e 1 e 2 e 2 


t _ t 


a 2 f 'f 2 <v(U »(!;)> l + (12) 

CD JJ^ \ 2 12 


0 0 


In view of Equation (II), take t = t and t = t + T . Equation (12) becomes: 

1 2 


<x*>=e- 2 ^e-^Va’/df / +T d 5 R « -«> e 0 ' <V + V (13) 

ro vO 1 «/o 2 M 2 l 


where R ( 5 - I ) = ( M (5 ) H (5 )^ (recal I that ( (J 2 > = 1 ). The form of Equation (1 3) 
H 2 1 1 2 

suggests the obvious change of variable to center of mass and relative coordinates: 


5, + * 

1 2 

2 


I = ? - 5 

2 1 


3 ( 5 , n 

The Jacobian J = ^ ,J jnr 

3 V 17 


1 so that d£ d£ = d£ d£ 

0 1 2 


The range of £ is from -t to t + T . If the £ integration is performed first, this range 

0 

gives the integration limits for £ , and the limits for £ become a function of £ as well as 

0 

t and T . The limits on £ are seen with the aid of Figure 3, a sketch of the integration 

/ \ 0 


domain in|| , £ |space. These limits 


5 < T , 


if < f o s i + 4 


I > r , 


-UP- < - < t +t - -i 
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Equation (13) becomes: 



where 


I 

1 


I 

2 


/ 


t + I 

2 0 " 2 ^ At 

e d£ 


LU 


2,3 t 

~w 


e^-e^l e“ 2 ^ 


(18) 


1 


t + T 

2 At 

e o a£ 


ILL 


23t 

k 

~28~ 


2 3r -3| 3|l| -23 1 

e e - e 1 1 e 


Using I and I in Equation (17) and re-arranging terms, 
1 2 
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(xx > = 


€ 2 a 2 
ao 


23 


l /: 




t + r 


~ K R © - 

r 


— e 


-2/3t 


/ 


t + r 


R (?) d? 
F 


(19) 


00 

Now, R (?) — 0 for ? » 7 = /* R (?) d? . (Note that ? must simply be larger than 

H M y Q J 

several T for this condition to be satisfied.) If we assume that R — ► 0 quickly enough so 
at F I 1 

tl._i n / *»\ r* .1 f i ~w- I l* • i_ j. J j. i -r* L.. 


~ p 

that e M ^ R (?) — 0 , then for t » T we may replace the limits -t and t + T by - 
H H 

and +oo • Doing so, and taking the limit t— oo so that the last term in Equation (19) 

• i 


< X x> 


< 2 o . 2 

00 -3t 

~w~ e 


R m (?) d? + e 2(Sr f e“^ R^?) d? 

-oo T 


/' 


( 20 ) 


If 1 /j3 » , then e^ and e ^ may be replaced by 1 in the integrals. Making this 

approximation, and re-arranging the integrals, we obtain: 


< * *> - 4^ I (* 3r+ ^ ) f y«> T )f y «> 


d? (21) 


For T y> t , the upper limit on the second integral may be changed to oo . It is expected 
F 

that 1/(3 » T (and this approximation has been made above), so that t will be larger than 

F , ^ 

r over most of the range of the correlation function. This gives / R (?) d? = t for 

H Jo F s s F 

both integrals, so that the result for T » r (again, this condition is just r larger than 

F 

several t ) is 
F 
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2 _ 2 


< x x > = e a 


M ~@ T 


oo /3 


The correlation function is: 


R (T) = = e- ffT 

* <x 2 > 


The integral scale of x fluctuations is thus 


/ 


oo 


-0T _ 1 

e dr - - 


( 22 ) 


(23) 


(24) 


This bears out the intuitive notion that ,3 is closely related to the integral time scale of x 
fluctuations. 

The mean square x excursion is given by using this in Equation (8): 


< x 2 > = £ 2 a 2 t t 

N ' co H x 


(25) 


For T< T , Equations (22) through (25) are not valid. For T<T , (3 t <5C 1 , so that 
M M 

-f* 

e~ & 1 + | 3t . Equation (21 ) may be written: 


<xx> 


e 2 a 2 

oo 


2/3 


oo r ) 

J (?) d? ' 2 Brf R m ({) df 


(26) 


For T « T , the second integral in Equation (26) is approximately T. Thus 
H 
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( x x ) r: 


1 


(27) 


2 2 

€ a 

CD 


/5 


for small t. This result will be used later in discussing the predicted spectrum. 

2.1 .3 Pressure Flucutations — If the pressure due to the shock wave at its mean location is 
P(x) , then the pressure at time t is given by P (x - x(t)). The fluctuating pressure from 
shock oscillation is then 


P (x , t) = P (x - x (t) ) - P (x ) 


(28) 


where x denotes the streamwise coordinate and x(t) is the excursion. The mean square 
pressure fluctuation is: 

< P 2 > = < [ P (x - x(t)) - P(x)] 2 > (29) 


If £ Is small, P (x - x) may be expanded In a Taylor series: 

<P 2 > = < [P(x)-x(t) P'(x) + .... -P(*)] 2 > (30) 

p 

so that for £ <K -^ 7 - 

< P 2 > = <[l (t) P' (x) ] 2 > = [ P' (x)] 2 <x 2 > (31) 

1 x 

This is accurate only when <( x 2 2 is small compared to shock thickness. If < x 2 2 is 

not small. Equation (29) would have to be used. This would require the probability distribution 
function of x , which has not been calculated. 
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2.1 .4 Spectral Density — The spectral density of the mean square pressure fluctuations is 
defined as the Fourier transform of the correlation function: 

r 00 

0(f) = 4 < p 2 > / R (t) cos 27rf dr (32) 

J o p 

where the factor 4<p 2 > appears so that Z* 00 <f> (f) df = < p 2 y . This definition of 0 is 

Jo 

consistent with the notation of Reference 1 . 

Assuming that < x 2 > is small enough so that Equation (31 ) is valid, R (t) = R (t) . Using 

P x 

Equation (23) in Equation (32), it is straightforward to obtain 


0(f) 


0 


4 <P 2 > 



2 


(33) 


Because Equation (23) is not accurate for r< T , Equation (33) is not accurate for 

P 

f> ~27TT~ * For f » 2tt~ ' ' s 9' ven ^ Y transform of EquaHon (27). This 

nr v nr v 

behaves asympt-otical ly as f 3 , compared fo f 2 for EquaHon (33). EquaHon (33), then 

will tend to be too large at high frequencies f > . 

Z7TT 

M 

2.2 Prediction of Pressure Fluctuations Ahead of a 45° Wedge 

The fluctuating pressure intensity and spectra are given by Equations (31) and (33), with 

Equation (8) giving < x 2 > and EquaHon (24) relating T t o 0 . Given the mean flow 

conditions, the parameters required are the boundary layer turbulent intensity e and integral 

scale t , and the integral scale r of the shock motion. 

H x 

Predictions will now be made for the 8 inch, 45° wedge of Reference 2. This particular set 
of measurements is chosen because the two-dimensional geometry is simple, and because of the 
good quality of both fluctuating and mean measurements. 
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2.2.1 Flow Parameters — The following are the mean flow parameters required: 


M = 2.0 

CD 

6 = 4 inches 

o 

h =8 inches 

The pressure gradient at the mean shock location is required. Figure 4 gives the mean pressure 
coefficient and root mean square fluctuating pressure for several test models in Reference 2. 
Measuring the slope for the 45° wedge case at the mean shock location, 

q q 

P'(0) = 0.42 - 2 - = 0.21 (34) 

o 


From Equation (32), the integral scale is given by: 


T. 

i 



0.(0) 

v> 


(35) 


where i is p or x. The integral scales are easily found from the 0 frequency point in 
Figure 2, and are 


6 

t = 25 — i 

x u 

0 


6 

r = 1 .9 — - 

H u 

0 


(36) 


i 

From Figure 4, < p 2 > 2 / q^ = 0.003 for the boundary layer. Turbulent velocity fluctuations 

are related to pressure fluctuations by (Reference 6): 


<p 2 >* = 0.7 p <u 2 > = 0.7 p € 2 a 2 (37) 

00 
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Taking p = p , 
00 


M 2 

oo 

1 .4 


<p!> 


oo 


(38) 


i 

Using < p 2 > 2 = 0.003 q , 

oo 


e 2 = 0.00214 M 2 (39) 

oo 


All required parameters have now been obtained. 

2.2.2 Root Mean Square Pressure Fluctuations — Using the above values in Equation (8), 


<x 2 > 


2 2 

€ a 

00 


r 


H 


T 

X 


= 0.00214 


6 

M 2 a 2 . 25 
00 00 u 

0 


• 1.9 


8 


o 


u 

0 


= 0.102 8 2 
o 


so that the root - mean square displacement is 


<x 2 > 2 = 0.32 8 (40) 

0 

This is smaller than the shock thickness (which is of order 8 ), so that Equation (31 ) may be 

used for the pressure. The fluctuation at the mean shock location is 
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( 41 ) 


<p 2 > 2 = p. (0) <x 2 > 5 = 0.214 -0.32 8 

o 


<P 2 > ? 


= 0.68 


oo 


This is in excellent agreement with the value of 0.7 in Figure 4. 


2.2.3 Spectrum — Figure 5 shows the predicted spectrum, Equation (33), along with the 
shock spectrum from Figure 2 and two spectra measured by Robertson (Reference 1 ) for three- 
dimensional protuberances. Except for being slightly high at higher frequencies, the agreement 
is excellent. fS/u = 1 0 -1 corresponds approximately to f = ' ' w ^ ere Equation 


H 

(33) is not expected to be accurate. The slope of the measured spectra at this point corresponds 
to f 2,6 , which is between the f -2 behavior of Equation (33) and the f -3 behavior of the 
high frequency limit discussed in Section 2.1 .4. A spectrum calculated from Equation (21) 
would show even better agreement with the data. 


Above fS /u w 2 * 1 0 -1 , the measured spectrum is dominated by separated flow 
u oo 

fluctuations, as pointed out by Robertson (Reference 1 ), and the present analysis does not 
apply. 
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Figure 5. Comparison of Predicted Spectrum with Measured Spectra 
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3.0 CONCLUSIONS 


A simple mechanism has been proposed for the pressure fluctuations near a shock wave inter- 
acting with a turbulent boundary layer. The shock wave is convected from its mean location 
by velocity fluctuations in the turbulent boundary layer, while stability of the mean flow tends 
to restore the shock wave to its original position. The restoring mechanism is assumed to be 
linear. The motion of the shock wave, given by Equation (5), is quite similar to the Brownian 
motion of a particle. The mean square displacement tends to a constant value at large time, 
given by Equation (8), so that the motion may be treated as a stationary random function. The 
mean square pressure fluctuation. Equation (31), and spectral density. Equation (33), calculated 
from this shock motion are in excellent agreement with experimental data. 

Because of the excellent agreement of the present theory with experimental data, and the 
straightforward physical model employed, it is concluded that turbulent boundary layer 
fluctuations are the dominant cause of shock wave oscillation in the case of flare and pro- 
tuberance induced separated flow. 
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